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Electrons hopping on the sites of a two-dimensional Lieb lattice and three-dimensional edge cen- 
tered cubic (perovskite) lattice are shown to form topologically non-trivial insulating phases when 
spin-orbit coupling is introduced. These simple models on lattices with cubic symmetry show a 
Dirac-like structure in the excitation spectrum but with the unusual feature that there is a dis- 
persionless band through the center of the spectrum and only a single Dirac cone per Brillouin 
zone. 



I. INTRODUCTION 

The study of topological insulatorJ^'^l ^as taken off in 
the last few years, proving to be an exciting area of 
research in condensed matter physics, while also hav- 
ing ties to fundamental physics due to the possibility 
of providing a testbed for exotic particles like Maiorana 
fermions,'^ fractionally charged vortices,!^ axion^SE and 
magnetic monopoles.l^ 

In general, topological insulators are bulk insulators 
having an energy gap between the valence and conduc- 
tion bands, but with gapless edge states (2D) or surface 
states (3D), caused by spin-orbit coupling, which are im- 
mune to non-magnetic impurities and geometric pertur- 
bations. Similar in spirit to the quantum Hall effect, 
the novelty of a topological insulator (TI) lies in the fact 
that it can be characterized by a topological invariant, 
and is not the result of a spontaneously broken symme- 
try. However, the invariant can not take on any integer 
value, as in the in the quantum Hall case, but is instead 
a so called Z2 topological quantum number that we label 
V, which can be either or 1. In 2D a single Z2 invari- 
ant is sufficient for the job, whereas in 3D one needs 4 
Z2 invariants {i'o',viV2i^3), which distinguish 16 possible 
topological phases and fall into the two general classes: 
weak (WTI) and strong (STI) topological insulators.^ In 
3D if there are an odd number of surface states then we 
have I'D = 1 and the system is in the STI phase. If there 
is an even number of surface states (possibly zero) then 

— Q and wc have a trivial insulator or WTI if any of 
the Vi are nonzero. 

The 2D topological insulator, also known as the quan- 
tum spin Hall (QSH) phase, was independently proposed 
by two groups and to date a variety of models have 
been established that can support the TI phase. The 
original model on the honeycomb lattice was realized 
by superposing two copies of the Haldane modelj^ one 
for spin up and other for spin down electrons, with the 
spin up and down electrons moving in opposite directions 
along the edge.l^ It was then predicted theoretically,'^! 
and later confirmed experimentally,^-' that this effect 
is present in the HgTe/CdTe quantum wells. Fu and 
Kane extended the concept to 3D with a toy model on 
the diamond lattice,^ and also predicted that the alloy 
Bii^j^Sb^:, which has large spin orbit coupling, would 
be a three dimensional TI.'^ This was later verified ex- 



perimentally by ARPES measurements which directly 
observed the topological surface states.fi^ Other models 
with a low energy Dirac structure, namely the 2D kagome 
and 3D pyrochlore systems, have been recently shown to 
support topological states.'^^H^Sl Experimentally the most 
promising system is currently Bi2Se3, which is known to 
have a large bandgap and a single surface Dirac cone.l^. 
Although not experimentally verified, there is also some 
theoretical work on to polog ical insulators in cold atoms, 
which looks promising.l^SM] 

On a practical level, these gapless surface states are 
robust against disorder, since time reversal symmetry 
disallows back scattering from non magnetic impurities. 
When gapped by time-reversal symmetry breaking per- 
turbations, superconducting or excitonic pairing, these 
surface states give rise to interesting insulating phases 
with exotic quasiparticle excitations already mentioned 
above. Electric manipulation of the spin degrees of free- 
dom in topological insulators should also be possible and 
they can therefore be of value in spintronics applications. 
Some of the existing topological insulators are also known 
to be among the best thermoelectric materials, such as 
the material Bi2Te3 discussed by Zhang et. dl.^^ and 
thus discovery of new materials or further understanding 
of existing ones would be of value. 

In this paper, we present another possible host for a 
topological insulator, namely the Lieb lattice in 2D and 
its 3D counterpart the perovskite or edge centered cu- 
bic lattice. We show that a simple tight-binding model 
for electrons on these lattices leads to TI behavior in the 
presence of spin-orbit coupling. These toy models are 
novel in that unlike most of the existing models, they 
have a simple cubic symmetry, and exhibit only a sin- 
gle Dirac node in the spectrum intersecting a degenerate 
fiat band precisely at the Dirac point. We also briefiy dis- 
cuss other perturbations that lead to topologically trivial 
gapped phases. 



II. LIEB LATTICE 

We begin with the tight-binding Hamiltonian 

Ho = -t ^ cLc,„ -I- h.c, (1) 
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FIG. 1: (a) Lieb lattice showing 3-site basis in unit cell, NNN 
hopping (red dotted lines) and basis vectors ai and a2. (b) 
Tight-binding dispersion for "H^- (c) Bandstructure for strip 
of Ny = 10 unit cells with open boundary conditions along 
the y direction and infinite along x with a — 0.2t. (d) As in 
(c), but with A = 0.2t instead. 



where creates an electron of spin a on site i of the so 
called Lieb lattice, seen in Fig. [T^, and t is the hopping 
amplitude for nearest neighbour sites. 

In momentum space Eq. ([I]) becomes Ha = 

J2ka *L'^k*k<^ with *ko- = (cik(T,C2k<T,C3k<T)^ and 



nl = -2t 



cos(A:2;) cos(fcy) 





where the lower triangle of the matrix is understood to 
be filled so that the matrix is Hermitian. The spectrum 
of "Hk consists of one degenerate flat band e'^'' = and 
two dispersive bands 



(1.2) 



= ±2t^/ cos^ kx + cos^ k 



(2) 



where the Brillouin zone spans — f<fca;<f,— f< 
fey < f • Plotting these bands yields a low energy Dirac 
like dispersion (seen in Fig. [Tja) , with a single Dirac cone 
in the Brillouin zone. At one third filling then, band 1 will 
be completely filled, with the degenerate flat band and 
the upper band being empty, and this state will behave 
as a gapless band insulator. 

There are several perturbations to the basic Hamilto- 
man ^ that open up a gap while respecting the trans- 
lational symmetry of the lattice. The simplest option is 
to add an on site energy e to all the sites that have only 
two nearest neighbours in Fig. [l^ above, or conversely 
the inequivalent sites with four nearest neighbours. We 
note that this set up has been discussed recently in the 
context of ultracold fermionic atoms. Here, one of the 



dispersing bands remains touching the flat band while the 
other dispersing band becomes isolated. Another option 
is to add a dimerization term, which staggers the hop- 
ping amplitude along both the x and y directions so that 
t — >■ t ± a. In this case the flat band becomes isolated 
from the dispersing bands, which are gapped symmet- 
rically (shown for strip geometry in Fig. [ij:) . Lastly, 
including a Rashba spin orbit term will again isolate the 
flat band, while also splitting the spin degeneracy among 
the valence and conduction bands. In all of the above 
cases the resulting insulating phases are characterized 
by conventional broken symmetries and are topologically 
trivial. 

The recent paper by Green, Santos and ChamorPsl 
discusses a similar set of energy bands having a flat 
band directly through the middle of a single Dirac cone, 
which have been created via staggered flux phases on the 
kagome lattice. Although their interests lie with the flat 
band itself, they discuss possible insulating phases that 
can be introduced in order to isolate the flat band and 
determine that, for spinless fermions, time reversal sym- 
metry must be broken in order to achieve this. Seeing as 
the Lieb lattice requires no magnetic field to create the 
same energy band structure, perhaps it is of some interest 
then that we can isolate the flat band without breaking 
time reversal symmetry by way of the dimerization term 
mentioned above. 

Now, to see if this model can support a TI phase we 
include the intrinsic SO interaction term 



so 



iX 



E 

{{i3))aP 



(3) 



where A is the amplitude for the next-nearest neighbour 
spin-orbit-induced interaction (shown in Fig. [l^), the 
term i/jj = (d-j x dfj)2=±l, where d\j and d?, are the 
two unit vectors along the nearest neighbour bonds con- 
necting site i to its next-nearest neighbour j and cr is 
the vector of Pauli spin matrices. This will also lead to 
the formation of a gap at the Dirac point while preserv- 
ing T and the translational symmetry of H^. Fourier 
transforming, we have 



= ±4A 





-i sin (fcj.) sin {ky ) 




(4) 



where the +(— ) sign refers to spin up (down) elec- 
trons, and the spectrum for the full Hamiltonian, Hk = 
then consists of the same degenerate flat band 



bands 



n/SO 



0, and the modified doubly degenerate dispersive 



E, 



(1:2) 



±2 



4A2 



n.y, 

(5) 

with a gap Ago = 4|A| at the Dirac point. 

We note that the flat band eliminates the possibility of 
writing down a low energy Dirac expansion for the system 
in terms of the Pauli matrices, as discussed e.g. in Ref. 
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[TSland other topological insulator models, but, following 
Ref. [26], it is possible to express the full Hamiltonian 
above in terms of a set of 3 x 3 matrices that form a 
spin-1 representation of SU{2). Although not a focus 
here, this makes it possible, for instance, to calculate the 
Chern numbers of the bands analytically. 

To prove that our system is indeed a topological in- 
sulator, we now show by an explicit calculation that it 
possesses a nontrivial Z2 invariant. There are several 
ways to do this in practice, although they are found to 
be equivalent in the end, so we use the method that is 
most convenient. 

According to Fu and Kan^^i when a crystal possesses 
inversion symmetry the Z2 topological invariant v is re- 
lated to the parity eigenvalues ^2m{^i) of the 2m-th oc- 
cupied energy band at the four T-invariant momenta 
r^. Our system is inversion symmetric and so we can 
use this method to find f. If we select site 1 of the 
unit cell as the center of inversion then the parity op- 
erator acts as 'P[il!i{r),tp2{r),ip3{r)] = [V'i(-r), ?/^2(-r - 
3-1 ): V'3(— r — 8-2)] on the triad of the electron wavefunc- 
tions in the unit cell labeled by vector r. In momentum 
space the parity operator becomes a diagonal 3x3 matrix 
= diag(l, e"*'*!''', e"*'^^''') and the four T-invariant 
momenta can be expressed as Ti — T:[xni + ymi)/2 with 
ni,mi — 0, 1. The eigenstates of can be found an- 
alytically in this case and it is then easy to determine 
the parity eigenvalues of the occupied bands. At | filling 
we find that three ^'s are positive and one is negative. 
Which of the four ^'s is negative depends on the choice 
of the inversion center but the product \liS,{Ti) = (—1)" 
is independent of this choice and determines the non- 
trivial Z2 invariant 1/ = 1, confirming that the system is 
indeed a topological insulator. Similar considerations for 
I filling also yield v = 1. 

Next, we solve the model in a strip geometry numeri- 
cally using exact diagonalization, in order to incorporate 
edge effects and demonstrate the bulk-boundary corre- 
spondence, which states that whenever v = 1 there will 
be a pair of topologically protected gapless modes along 
each edge in the system. For the trivial insulating state 
with a = O.lt, we find no edge states (Fig. [l]:), but for 
A = O.lt, we indeed find a pair of spin- filtered gapless 
states associated with each edge, as seen in Fig. [TJl. As 
there are 3 atoms in the unit cell, the two edges of the 
slab are not equivalent, and therefore the edge states are 
not degenerate. 



III. PEROVSKITE LATTICE 

The perovskite lattice, depicted in Fig. [2^, is a 
straightforward generalization of the Lieb lattice into 3D 
and can be viewed as a simple cubic lattice with addi- 
tional sites positioned at the centers of all edges. 

As with the Lieb lattice, our starting point is the tight 
binding model given in Eq. ([T]), which in momentum 
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FIG. 2; (a) Perovskite lattice showing 4 sites in unit cell 
along with basis vectors, (b) High symmetry points in the 
Brillouin zone, (c) Bandstructure inside bulk along path of 
high symmetry for (d) As in (c), but with a = Q.2t. 



space gives 



nl = -2t 



I cos kx cos ky cos k^ 


V 



The spectrum then consists of two degenerate fiat bands 



(3,4) 



and the two dispersive bands 

(1,2) 



±2iA/cos^ fcj, -t- cos^ fcj^ + cos^ /cz. (6) 



Again, we have a single Dirac point in the Brillouin zone 
and at \ filling the system behaves as a gapless band 
insulator. The spectrum is shown in Fig. [2]:. As above, 
we can create a topologically trivial insulating phases by 
including an anisotropic on site term or a dimerization 
(Fig.[2}i). 

To open up a topologically non-trivial gap, we con- 
sider spin-orbit coupling Eq. pk. Since the d^j^ now lie 
in three-dimensional space, the~Hamiltonian does not de- 
couple for the two spin projections as above, and instead 
becomes an 8 x 8 matrix in /c-space. The spectrum can 
be found numerically and is shown in Fig. [3^. Unlike the 
model on the Lieb lattice, the addition of the spin orbit 
term creates a dispersion in the degenerate flat bands and 
a gap that depends on the sign of A, which is no longer 
symmetric around zero energy. Switching the sign of A 
inverts the band structure. 

We now study the topological classes of these insu- 
lating phases. As above, the Zi topological invariants 
(y^;v\V2V-^ are easy to evaluate when a crystal pos- 
sesses inversion symmetry and can be determined from 
knowledge of the parity eigenvalues ^2m^i) of the 2m,- 
th occupied energy band at the 8 T-invariant momenta 
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FIG. 3: (a) Bulk bandstructure shown along path connecting 
high symmetry points for A = 0.2t. (b) Bandstructure for 
slab geometry along lines connecting four surface TRIM for 
A = 0.2f. 



(TRIM) Ti that satisfy = + G. The 8 TRIM in 
our system can be expressed in terms of primitive recip- 
rocal lattice vectors as ^ i={nxn.in2,) = ("ibi + «2b2 + 
n3b3)/2, with rij = 0, 1. Then is determined by the 
product (-1)''" = n„,=o,i'5nin2«3. and (-l)--i^^.3 = 

n„,^,=o,i;«.=i '^"i«2n3, where Si = 11™=! 6™ (r^)- 

If we select site 1 of the unit cell, Fig. [2] as 
the center of inversion then the parity operator acts 
as 7'[V'i(r), V2(r), V'3(r),V'4(r)] = [V'i(-r), ■02(-r - 
Si):V'3(~r ~ a-2),'04(— r — as)] on the four-component 
electron wave function in the unit cell labeled by vec- 
tor r. In momentum space and including spin the par- 
ity operator becomes a diagonal 8x8 matrix — 
diag(l, e'"^!-'', e-*'^^ '', e~''^3 '')(X)diag(l, 1). It is straight- 
forward to obtain the eigenstates of "Hpi and the parity 
eigenvalues of the occupied bands numerically, then de- 
termine the Z2 invariants. At quarter filling, we find that 
(5 = -1 at the iJ point and 5 = 1 at the other TRIM, 
so the spin-orbit phase is a (1; 111) strong topological 
insulator. 

Incorporating edge effects into the system with a slab 
geometry and plotting the band energies along lines con- 
necting the four surface TRIM (Fig. [3]d), we can also 
see the bulk energy bands and an odd number of surface 
states which traverse the gap, a behavior characteristic 
of STL 



IV. CONCLUSIONS 

We have established a model for a topological insula- 
tor on lattices with simple cubic symmetry. In 2D we 
have shown that a tight-binding model with spin-orbit 
coupling on the Lieb lattice supports gapless topologi- 
cally protected edge modes, and possesses a non-trivial 
Z2 invariant. In 3D, we have demonstrated that a similar 
model on the perovskite lattice supports 2D gapless sur- 
face states when in the (1;111) STI phase, brought on by 
spin orbit coupling. We have also identified other gapped 



phases in these models that are topologically trivial. 

In addition to the topologically non-trivial band struc- 
tures, the models considered here exhibit flat bands, 
which could give rise to interesting strongly corre- 
lated electron states even for relatively weak interaction 
strength. 1221 A necessary condition for non-trivial corre- 
lated physics is that the single-electron states forming 
the fiat band are delocalized in space (localized states 
typically lead to formation of a Wigner crystal). Delo- 
calized single-electron states are guaranteed to exist if 
the fiat band possesses a non-zero topological invariant, 
as happens e.g. in quantum Hall liquids. Unfortunately, 
this is not the case in models considered in this study. 
In the Lieb lattice spin-orbit coupling preserves the flat 
band and separates from other bands. The resulting flat 
band however turns out to be topologically trivial, i.e. it 
has = 0, while the TI behavior derives from the two 
dispersive bands which have v — \. In the perovskite lat- 
tice spin-orbit coupling produces a significant dispersion 
in the original fiat band and this will limit the importance 
of correlations in the system. 

Can our model be realized in a physical system? There 
exist many perovskites in nature as well as many layered 
perovskites composed of weakly coupled 2D planes with 
Lieb lattice structure. The most prominent example of 
the latter are the Cu02 planes in high-Tc cuprate super- 
conductors such as YBa2Cu307 or Bi2Sr2CaCu208. In 
3D, our model system is an idealization of the naturally 
occurring perovskites, most of which also have a heavy 
central atom in the middle of each cubic cell, such as 
SrTiOa or the double perovskite structure Ba2NaOs06. 
In these real materials the electron behavior near the 
Fermi level derives from the eg and tig orbitals of the 
transition metal element occupying the cubic site while 
the edge sites are normally oxygens whose p-orbitals are 
far away from the Fermi level. The resulting band struc- 
ture for the active orbitals is then significantly more 
complex than that captured in our simple tight binding 
model. Nevertheless, our model calculations demonstrate 
that this class of tight-binding Hamiltonians on lattices 
with cubic symmetry can support topological phases, 
both in 2D and 3D. We hope that our work will stim- 
ulate detailed band structure calculations of perovskites 
with heavy transition metal elements in search for new 
families of topological insulators. 

In the more exotic realm, it might be possible to ar- 
tificially engineer the 2D system by modulating a two- 
dimensional electron gas with a periodic potential having 
Lieb symmetry, as achieved recently in constructing the 
'artificial graphene'.'^ Another possibility lies with cold 
Fermionic atoms in optical lattices as discussed in Refs. 

Eng. 
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